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A simple electrostatic theory of condensation of rod-like polyelectrolytes under influence of poly- 
valent ions is proposed. It is based on the idea that Manning condensation of ions results in formation 
of the Wigner crystal on a background of a bundle of rods. It is shown that, depending on a single 
dimensionless parameter, this can be the densely packed three-dimensional Wigner crystal or the 
two-dimensional crystal on the rod surfaces. For DNA the location of charge on the spiral results in 
a model of the one-dimensional Wigner crystal. It is also argued that the Wigner crystal idea can be 
applied to self-assembly of other polyelectrolytes, for example, colloids and DNA-lipid complexes. 
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Many -rod-like polyelectrolytes, such as double he- 
lix DNAEJ, F-actin, microtubules, and tobacco mosaic 
virusa are known to self-assemble into bundles of parallel 
densely packed rods. All these macromolecules are neg- 
atively charged and their lateral association is induced 
by Z-valent cations where Z > 2. For DNA this phe- 
nomenon is called condensation and is studied in con- 
nection with the dense packing of viral DNA. According 
to the mean field Poisson-Boltzmann theory, two paral- 
lel rods should always repel each other. Two physical 
phenomena which are not included in this theory were 

Suggested as possible reasons for the puzzling attraction 
liiJ. When distance between rods is large their attrac- 
tion is related to the correlation of thermal fluctuations of 
screening atmospheres of two rods. At smaller distances 
one should take into account that the Manning condensa- 
tion of ions on the surface of rods leads to strong spatial 
correlations of ions or even to their crystallization. In this 
situation, two crystals with proper phases attract each 
other. Below I calculate the binding energy of a dense 
bundle and therefore talk only about the second mecha- 
nism. I also resolve an important contradiction in the lit- 
erature. Most of the publications deal with two pods and 
calculate a pairwise force acting between j&emETEl. On 
the other hand, two recent publicationst30 claim that 
bundles are formed by non-pairwise-additive interaction. 

The first goal of this paper is to present a theory of 
attractive interactions of rod-like polyelectrolytes, based 
solely on electrostatic interactions. It considers forma- 
tion of densely packed bundles in a very dilute solution 
of cylindrical rod-like molecules of the radius r and length 
L ^> r. It assumes that the cylindrical surface of a rod 
is negatively charged with the charge density — e/b per 
unit length of the rod. Point-like positive ions with the 
charge Ze are added to the solution. 

The main idea of this paper is that a bundle of paral- 
lel rods can be considered as an uniform negative back- 
ground at which condensed ions form the Wigner crystal. 
The cohesive energy of this crystal is the reason for the 
attraction and the bundle formation. (This idea is similar 
to the theory of matter in a superstrong magnetic field 



of a neutron star where electrons behave as negative rods 
and nuclei form the Wigner crystal on the background of 
electron bundlea!3.) For the case in which ions can not 
penetrate the rods, the binding energy of the bundle per 
rod is calculated below as the function of the dimension- 
less parameter Zb/r. The result is that for Zb/r ^> 1, 
the three-dimensional densely-packed Wigner crystal is 
formed and interaction is not-pairwise-additive. In the 
opposite case, Zb/r <g; I the two-dimesional Wigner crys- 
tal appears on the surface of each rod and the interaction 
becomes short-range pairwise. My results in the first case 
qualitatively resemble Ref. II], while in the second case 
this paper is close to Ref. [Twhere idea of the surface 
Wigner crystal was originally suggested. None of these 
papers, however, produced simple analytical dependen- 
cies similar to those which I derive below using Zb/r as 
a large or small parameter of the theory. 

The second half of this paper concentrates on the 
specifics of DNA where surface charges form a spiral. 
I show that this leads to the appearence of the one- 
dimesional spiral-like Wigner crystal of ions on the sur- 
face of each rod. Interaction of such crystals determines 
the DNA condensation. It is also argued that the idea of 
a Wigner crystal can be applied not only to the rods, but 
also to the self-assembly of many other different polyelec- 
trolytes. I apply this theory to colloids and complexes of 
DNA with cationic lipid membranes. 

Returning to the main problem of rods and point-like 
counterions, I argue that for sufficiently large dimension- 
less parameter Z^, where £ = e 2 /bKksT and k is the ef- 
fective dielectric constant of the water-polymer system, 
rods condense in the maximum density cylindrical bun- 
dles of N parallel molecules. At large Z£ the charge 
of the bundle is almost completely compensated by the 
opposite charge of NL/b of positive ions which experi- 
ence the Manning condensation inside the bundle. To 
find a configuration of condensed ions, one can view the 
bundle as the uniformly charged cylindrical background. 
The potential energy of the Coulomb repulsion of ions in 
these conditions is much larger than their kinetic energy. 
Therefore, they form the Wigner crystal on the negative 
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background. Assume first that N is so large that the ra- 
dius of the bundle rN 1 ^ 2 is much larger than the lattice 
constant of the Wigner crystal R. In this case ions form 
a densely packed three-dimensional crystal. One can eas- 
ily calculate R from the condition that the charge of the 
Wigner-Seitz cell, R 3 /br 2 , equals Z. This gives 



R ~ (Zh 



.2x1/3 
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Comparing R with the bundle radius rN 1 / 2 one finds 
that the three-dimensional case takes place at N ^> N c = 
(Zb/r) 2 / 3 . (In the opposite case N <C N c ions form one- 
dimensional crystal as shown below). I assume that ions 
can not penetrate rods. Then ions can fit in the free 
space between rods without strong deformation of the 
three-dimensional crystal only at R 3> r. According to 
Eq. (0) this inequality means that Zb/r ^> 1. Now we 
see how the important parameter Zb/r appears in the 
theory. Assuming, for the beginning, that rods are nar- 
row and weakly charged so that Zb/r ^ 1 one finds that 
N c > 1 so that both situations N > N c and N < N c 
have to be considered. 

The binding energy of a crystal per ion, £j, can be es- 
timated as the energy of interaction of an ion with its 
Wigner-Seitz cell 



Z 2 e 2 
Ei ~ kR' 

Substituting R from Eq. ([[]) one obtains 
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e .^^ Wr )2/3 (Zb/r ^ I; iV»A c ). (3) 
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If N <C N c , the bundle is narrow and the ions form a one- 
dimensional Wigner crystal, lattice constant of which, R, 
is equal to Zb/N . Then, the energy of interaction of an 
ion with its Wigner-Seitz cell is 



Ze 2 N 
Kb 



ln{N c /N) {Zb/r > 1; N < N c ). (4) 



It is easy to verify that Eq. (0) matches Eq. (|j) at 
N ~ N c . Thus, at Zb/r 3> 1, the energy per ion, 
Si, grows linearly at small N and reaches saturation at 
N 3> N c . This is a clear demonstration of a non-pairwise 
additive interaction. Of course, the surface correction to 
the energy of the bundle provides some growth even at 
N N c and leads to formation of macroscopic bundles. 

In order to obtain the binding energy per molecule, 
Eq. (|J) and Eq. (||) have to be multiplied by the number 
of ions per molecule, M = L/bZ. This transition does 
not change the dependence on N. Therefore below I will 
continue to present results in the form of £,. 

Before switching to the more complicated case Zb/r <C 
1, I discuss the condition on the temperature T or, in 
other words, on the parameter Zt; at which the suggested 
above theory is valid. Consider a large bundle of three- 
dimensional densely packed Wigner crystal, which stays 



neutral even at quite a small concentration of ions. If 
ksT <C £i the thermal motion can be neglected and my 
theory is valid. Using Eq. (||) one can rewrite inequality 
ksT <C £j as 



Ztl > (r/Zb) 
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At Zb/r 1 this inequality is much weaker than the 
standard condition Zt; 3> 1 of the Manning condensation 
of ions with charge Z at an isolated molecule (it follows 
from Eq. ([|) at N — 1), because of the simultaneous in- 
teraction of each ion with many rods. Even if one uses 
K = 81 of pure water, £ ~ 4 at room temperature and 
the condition of Eq. (g) is easily fulfilled. 

It is obvious that wnen Eq. (@) is valid, the maximum 
density bundle is more strongly bounded than one with 
the smaller density. Indeed, the decreasing density re- 
sults in the increase of the lattice constant R and, ac- 
cording to Eq. (|^), substantially diminishes £*. At the 
same time the increase of the entropy term in free energy 
can not compensate for this loss in the binding energy. 

As stated above, this theory works if fc^T -C £j. Ac- 
tually the Wigner crystal melts when ksT is yet numeri- 
cally much smaller than £,-, so in some experimental con- 
ditions or for low dimensionalities (see below) one deals 
with strongly correlated liquid, not a crystal. This, how- 
ever, does not change the estimate for its correlation en- 
ergy and for and leaves unchanged the validity criterion 
of the theory. The fact that it is the short range order 
of ions which is responsible for the attraction of rods at 
a finite temperature was emphasized in Ref. |3[ Below, I 
continue to use the Wigner crystal language because it 
creates a simpler image. 

Consider now more strongly charged and thicker rods 
for which Zb/r <C 1. In this case iV c <C 1 and, therefore, 
one has to deal only with a three-dimensional problem. 
However, the rods are so thick that their radius is larger 
than the lattice constant of the optimal densely packed 
crystal. If, as we assumed, ions can not penetrate rods, 
the optimal densely packed crystal can not be formed. 
Under this restriction, ions condense at the surface of 
each rod forming the two-dimensional Wigner crystalEI 
with the lattice constant 



R - (Zbr) 
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An insulated rod has a similar crystal at its surface. The 
binding energy of the bundle originates in narrow con- 
tact stripes where rods pairwise contact each other so 
that their Wigner crystals overlap. The width W of 
this stripe will be calculated below. Inside the stripe, 
densities of both negative background and positive ions 
are doubled and the local lattice constant R s becomes 
smaller (R s = R/y/2). As a result, according to Eq. (||), 
the cohesive energy of the crystal per ion of the stripe 
becomes larger than for the case of the two separated 
rods. Combining Eq. (||) and Eq. (^) one finds that the 
binding energy of the bundle per one ion of the stripe 
equals Z 3 / 2 e 2 /(6r) 1 / 2 . To get £j one must multiply this 



2 



energy by the fraction, p ~ W/r, of all ions of the bun- 
dle which reside in the contact stripes. The width W 
can be found from the condition that at the distance W, 
from the line of contact the surfaces of two contacting 
rods diverge from each other at the distance of the order 
of R. Indeed, the interaction between two-dimensional 
Wigner crystals is exponentially weak if the distance be- 
tween parallel planes in which they are situated is larger 
than their lattice constant R. A simple geometrical esti- 
mate gives W ~ (ri?) 1 / 2 . Finishing the calculation of £j 
one obtains 

7e 2 

£i~ —-{Zb/rf^ (Z6/r«l; AT » 1). (7) 
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Together Eq. (@) and Eq. (§ give binding energy per 
ion of a three-dimensional crystal at all values Zb/r. At 
Zb/r = 1 they obviously match each other. At Zb/r <C 1 
energy given by Eq. (m) is smaller than the one from 
Eq. (|3|). This is a natural result of the restriction that 
ions do not penetrate rods. 

Consider now the validity of the main assumption, that 
an array of negative discrete charges on the rod surface 
can be effectively replaced by an uniform negative back- 
ground. This idea works exactly only in the limit Z 1 . 
On the other hand, it fails at Z = 1 because in low tem- 
perature limit all ions and discrete negative charges form 
neutral Bjerum pairs and nothing depends on the mu- 
tual positions of the rods. Thus, £j = at Z = 1. What 
happens at Z > 2 depends on the spatial distribution of 
the discreet negative charges. If they are distributed ran- 
domly, then at Z — 2 Eqs. (||), (^) and ([?]) overestimate 
Si roughly by a factor of 2, while at Z > 3 they work well. 
On the other hand, when negative charges are clustered, 
Ei can be much lower even at Z = 3. Imagine, for exam- 
ple, that negative charges form compact triplets. Each 
one would be neutralized by an ion and gj would van- 
ish. More interesting effect of concentration of negative 
charges in DNA is discussed below. 

But first assume that the surface of DNA is uniformly 
negatively charged. For DNA b = 0.17nm and r = lnm, 
so that Z — 6 is the border between range of validity of 
Eq. (0) and Eq. ||. This means that for DNA in this 
model, attraction is due to the short range pairwise inter- 
action of molecules induced by correlation between their 
surface Wigner crystals. This conclusion is in agreement 
with Ref. |. 

Recall now that in DNA, negative charges are located 
along two spirals which are separated by wide and narrow 
grooves. For simplicity, assume that the width of the nar- 
row one is zero what leaves only one spiral of the negative 
charge with the double density. (This approximation, of 
course, overestimates the inhomogeniety of the negative 
charge distribution and can be easily avoided with minor 
change in the result.) Ions tend to concentrate on the 
same spiral and form the one-dimensional Wigner crys- 
tal along it. When two parallel DNA rods touch each 
other by these crystaline spirals (for this purpose one rod 
should be shifted along their direction by the half of the 



helical period) they create a spot with larger binding en- 
ergy per ion. Thus, attraction is produced in a way very 
similar to the case of contacting two-dimensional crystals 
discussed above. Extending this analogy for two contact- 
ing DNA molecules, one can introduce the contact stripe 
where attraction is generated. The only difference is that 
the width of this stripe W is smaller for DNA, because 
the Wigner crystal spirals of the two contacting DNA 
molecules diverge not only in the plain perpendicular to 
the rods but also along direction of the rods (spirals cross 
at a finite angle). One can find that W ~ Rr/nb ~ R, 
where R = nr/n is a period of the one-dimensional crys- 
tal and n — 10 is the number of base pairs in the helix 
period. Thus, roughly speaking, only one ion per helix 
period adjusts its position due to the contact and con- 
tributes to the binding energy of the two contacting rods. 
Multiplying Eq. (||) by the fraction of such ions R/2irr 
one obtains 



At Zb/r <C 1, Eq. (|J) gives smaller energy than Eq. (0). 
This is an effect of concentration of negative charges on 
the spiral. Note that this effect is exactly opposite to the 
prediction of Ref. [l^. The reason for this difference is 
that the authors of Ref. |l3| assumed that charge of ions 
is uniformly distributed over DNA surface, so that neg- 
ative and positive charges are strongly separated. Such 
a separation of charge costs a large electrostatic energy 
and can happen only if nonelectrostatic forces dominate. 
These forces are beyond the scope of my work. 

One might wonder whether each contacting pair of 
DNA molecules in a bundle can gain the energy given 
by Eq. (|^). This can be easily done in a square lattice 
of rods, with one square sublattice shifted along the rod 
direction by the half of the helical period. In a densely 
packed hexagonal lattice one can shift in the same way 
every second layer, so that each rod is attracted to 4 its 
neighbours. Thus in the approximation of one spiral both 
lattices have the same energy. 

Note that literally speaking this theory is applicable 
only to molecules of DNA which are shorter than its per- 
sistent length. It is known that in a weak DNA solution 
of longer molecules, each molecule self- assembles into a 
toroidal particle, where locally strands of all turns are 
parallel to each other. Theory of this paper is applicable 
to toroidal particles as well. Even longer DNA molecules 
condense in more complex globulesc-j. If such a globula 
consists of condensed bundles this theory can still work. 

Although three-dimensional densely packed Wigner 
crystal of ions is not realized in DNA bundles (for Z < 6), 
many other rod-like polyelectrolytes are known. Some 
of them are narrower or more weakly charged so that 
Zb/r ^> 1 and three-dimensional crystal of ions should 
appear. PPP3 with b/r ~ 2 is a good exampleEj. 

I would like to show now that the idea of Wigner crys- 
tal (or strongly correlated liquid) of ions for description 
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of their Manning condensation can be applied not only to 
rods but to quite a variety of systems. Imagine, for exam- 
ple, a rigid plane polyelecrolyte in the presence of counte- 
rions with a large charge Z. Such a system self-assembles 
in a dense stack of parallel sheets which provides uniform 
background for the Wigner crystal of ions. Once more, 
depending on the thickness of the sheets and their surface 
charge density, one can talk about the three-dimensional 
Wigner crystal of ions or about two-dimensional crystals 
in gaps between sheets. 

Counterions should not necessarily be small. The role 
of a countcrion can be played by a strongly charged rod- 
like polymer while the uniform background is provided by 
the stack of positive parallel sheets. Complexes of DNA 
with catipftic lipid membranes are good examples of such 
a systemlia. All the rods orient themselves parallel to 
each other and to the sheets. Their ordering in a plane 
perpendicular to their long axis depends on the sheet 
thickness. At a small enough sheet thickness, the circular 
cross-sections of rods are arranged in the two-dimensional 
Wigner crystal similar to vortex lattice in superconduc- 
tors, while at larger thicknesses they form parallel one 
dimensional crystals in the gaps between sheets. 

Returning to small counterions, consider the self- 
assembly of rigid spherical polyelectrolytes. Thinking 
about colloids, assume that the surface of a sphere is 
charged negatively and positive ions can not penetrate 
it. If spheres are large enough, ions form the two- 
dimensional Wigner crystal at the surface of a sphere. 
This leads to contact attraction of spheres similar to the 
one calculated above for two contacting cylinders. The 
binding energy of two spheres with the radius r and the 
charge —Qe, originates in the contact disc (an analog of 
the contact stripe for two cylinders) of the radius (ri?) 1 / 2 , 
where r is the sphere radius, R is the lattice constant 
of the Wigner crystal (see computer simulation of these 
discs in Ref. 0). The additional energy per ion of the 
disc is gained, because the disc has a denser background 
and a larger ion density than the rest of the surfaces of 
the spheres. Once more £j is given by Eq. (^). Multiply- 
ing this energy by the number of ions in such a disc, r/R, 
and substituting R ~ r(Z/Q) 1 / 2 one gets a surprisingly 
simple expression for the binding energy of two spheres 
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which looks like a Coulomb attraction energy of the 
spherical colloid particle and a single ion at its surface. 
Eq. (Js|) fails with the decrease of the ion concentration in 
the solution, c, when each sphere acquires the net nega- 
tive charge sufficient to overcome calculated attraction. 
This happens at fcsTZn(l/c) ~ Z 2 e 2 /nR. At lower tem- 
perature Eq. (^|) gives energy which is larger than pro- 
portional to T depletion energy of Ref. [l?]. 

In conclusion, this paper studies the attraction and 
self-assembly of rigid polyelectrolytes due to strong cor- 
relations in positions of polyvalent counterions. To de- 
scribe these correlations, I adopt the notion of the Wigner 



crystal. For rods I find that, depending on the value of a 
single dimensionless combination, Zb/r, optimal cofig- 
uration of ions on the background of a bundle is the 
three-dimensional densely packed Wigner crystal or the 
two-dimensional crystal on the surface of rods. When 
applying this theory to DNA condensation I introduce 
the notion of one-dimensional spiral Wigner crystal and 
estimate the binding energy of two DNA helixes due to 
the contact of these crystals. I also calculate the binding 
energy of two colloidal particles created by the contact 
of Wigner crystals of ions on their surfaces. Future de- 
velopments of this work should include such factors as 
behaviour of the dielectric constant of water at small dis- 
tances, the finite size of counterions, and more realistic 
distributions of polyelectrolyte charges. 
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